Isospin Squeezed States , Disoriented Chiral condensates and Pion Production: A 

Dynamic Group Theoretical Approach 
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I. INTRODUCTION 

Quantum optical analogies have long been use- 
ful tools in describing multipionproduction in 
hadronic collisions [J 0, IE 0, m- Quantities 
such as pion multiplicity distributions and corre- 
lations among identical and charged pions have 
been successfully explained on a phenomenolog- 
ical level by exploiting the similarity between 
light and bosonic particles such as pions. This 
has been very successful for systematizing data 
in high energy collisions since the "pre-QCD" 
days ,when the concept of coherent pion produc- 
tion was introduced p, 0, E[ . Since then, many 
distributions which have their origin in photon- 
counting experiments have found an application 
in particle physics p. Inparticular the nega- 
tive binomial distribution |2|, the Perina-McGill 
distribution kjL and the squeezed coherent dis- 
tribution 0, El have proved very useful in var- 
ious high energy applications. Recently these 
quantum optical states have been revived as pos- 
sible states for the disoriented chiral conden- 
sate, which is conjectured to occur in relativistic 
heavy ion collisions 0,00. 

The distinguishing factors between optical and 
pionic physics are the conservation laws and fi- 
nal state interactions between pions 0. Thus, 
in translating the quantum optical analogy to 
pion production, care has to be taken to pre- 
serve the conservation laws such as isospin, 
charge etc. obeyed by the underlying inter- 
actions 
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This is done by impos- 
ing constraints which exploit the group theo- 



retic structure of the symmetries and involve 
projection operator techniques [20]. In partic- 
ular for strong interactions where isospin is con- 
served, coherent states of pions are obtained by 
globally by projecting out a particular isospin 
state from a general number operator coherent 
state. Such a global projection assumes that all 
pions are assigned the same momentum func- 
tion and thus isospin dependence of event by 
event quantities like momentum dependence of 
pion correlations (HBT effect) and back-to-back 
particle-antiparticle correlations (PAC) cannot 
be described in this approach. One has to do a 
full quantum field theoretical treatment in mo- 
mentum space and find momentum dependent 
isospin variables to give a full event by event de- 
scription of the isospin coherent state. 

Recently a number of papers in literature have 
addressed the question of isospin conservation in 
the context of pion production associated with 
the disoriented chiral condensate (IE 21 , 2^, 3(3 . 
These involve the construction of the isospin 
squeezed states. However, these papers do not 
make any proposal for the dynamical origin for 
the anomalous production of particles in relation 
to the squeezed states. This has motivated us to 
look deeper into the question of whether isospin 
squeezed states can be dynamically generated in 
models for pion production in certain high en- 
ergy processes. Furthermore, only the isospin 
structure at zero relative momentum of pions 
has been investigated (2l|. Our approach is the 
most general one which has the results of these 
papers as special cases, and is more suitable for 



2 



the dynamical origin of pion distributions seen in 
physical processes. It can also be easily general- 
ized to include disoriented chiral condensates., so 
that isospin invariance in DCC production does 
not have to be compromised. 



II. THE PION HAMILTONIAN 

The foundations for the model used are laid in 
our previous paper j23( , hitherto referred to as 
I. In the next two sections , to make this paper 
complete and also to relieve the reader from the 
ardrous task of continually referring to I, we out- 
line the salient features and results derived in I. 
Since this is the basis for our subsequent discus- 
sion and results, the repetition is not redundant 
but useful. 

In I, we used a background field analysis of the 
0(4) sigma model with symmetry breaking keep- 
ing one-loop quantum corrections to construct 
the most general Hamiltonian for the evolution 
, formation, decay of the DCC in an external, 
expanding metric. 

The complete action for the pion-sigma part 
of the 0(4) sigma model which preserves chiral 
symmetry is given by 

S = J d 4 x^(d^d^+d^Z)~(7? 2 +Z 2 -fZ) 2 

(1) 

Where, the chiral field is 



linear in the sigma field. The action becomes 



$(r,t) = £(r,t) +it 
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r,t) 



(2) 



The true vacuum is < $ > 2 =< £ > 2 + < 
7r > 2 = f 2 and is taken traditionally as < E >= 
/■^and < 7i >= 0. The choice of this vacuum 
implies that chiral symmetry is spontaneously 
broken. Note that the fermionic part of the 
original sigma model Lagrangian has been ne- 
glected here. This is allowed since our focus is 
on the condensate formed in the symmetry bro- 
ken phase where quark degrees of freedom are 
already confined. The explicit breaking of sym- 
metry is implemented in the Lagrangian to give 
the pions mass through the presence of a term 



S 



d x 



v 2 ) 2 + e£ 



(3) 

where f n is replaced by a general parameter v , 
which in the limit e — > goes to /„-. Requiring 



the minimum to be still (0, /„. 
we have 



f 2 - 



A 



to leading order 



(4) 



and 



and the pion and sigma masses are m = 4- 

J TV 

m a 2 = 2\f 2 + m 2 n . 

f w = 92MeV is the pion decay constant and the 
meson masses are = 138MeV/c 2 and m a = 
600MeV/c 2 . 

Although , to give the pions mass, we have ex- 
plicitly broken the symmetry, if the explicit sym- 
metry breaking is very small compared to the rel- 
evant scale of QCD then it is still a good approxi- 
mation to apply the notion of spontaneously bro- 
ken symmetry, the symmetry breaking term tilts 
the double well potential given in (3). As long as 
the potential is tilted only slightly the pions are 
much softer than the sigma so the effect due to 
the spontaneous symmetry breakdown of chiral 
symmetry dominates the dynamics if e is small. 

The formation of the DCC takes place when, 
in a rapidly cooling and expanding plasma, the 
expansion rate of plasma is greater than the rate 
at which the field evolves from a state of restored 
symmetry to an equilibrium state of broken sym- 
metry 2 3]- This evolution of the field can be 
implemented by making the expectation value of 
$ a function of time. Such a study has been car- 
ried out in reference |23[ in great detail. An ad- 
ditional feature in the study of the formation of 
the DCC is that of the expansion of the plasma 
- two scenarios exist, the first being that of a 
sudden quench, where < $ > changes from to 
f n instantaneously and the second in which the 
system goes through a meta-stable, disordered 
vacuum typically represented by: 

E = UCos{6); TT = f n nSin(6). (5) 

The system then relaxes by quantum fluctua- 
tions to an equilibrium configuration. Here 9 
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measures the degree of disorientation of the con- 
densate. 

Some of the traditional signals of the DCC are 
the enhancement of low momentum pion modes 
llH and the anomalous charged to pion ratio of 
the soft pions [l^. The study of the DCC re- 
quires the further analysis of the pions produced 
by the decaying plasma and therefore it is imper- 
ative that the isospin structure , which is such 
an integral part of pion studies, has also to be 
incorporated. Studies carried out thus far have 
concentrated on classical isospin structure. We 
look at the isospin structure at the quantum level 
through the construction of an effective mean 
field Hamiltonian which allows us to examine the 
dynamical effects of isospin conservation. This 



provides us with a framework for examining the 
effects of isospin conservation on the pion multi- 
plicities. 

In |23| we have studied the effects of an so(4) 
sigma model with spontaneous symmetry break- 
ing in a spherically symmetric and homoge- 
neously expanding plasma. The line element is 
the FRW metric: 

d s 2 = dt 2 - a(t) 2 dx 2 , (6) 

where a(t) is the expansion parameter . We 
treated the quantum field fluctuation 
around the general classical background field pa- 
rameterized by three angles: 



/ fnCos(p) Sin{9) Sin(a) \ 
f n (p)Sin(9)Cos(a) 
f. K Sin(p)Sin(9) 
UCos{9) 



(v + \ 

V3 

V ° ) 



(7) 



where f+,f- are the vacuum expectation values 
of the charged pions and v% of the neutral pion 
. Without DCC formation these would be zero. 
In order to consider all the special cases that 
are possible in a transparent way, we simplified 
the parameterizations of the possible form for 
the background field to two angles, 9 and p by 
letting a = f . Then, v± = ^Cos(p)Sin{9), 

V3 = f n Sin(p)Sin(9), and a = f n Cos(9). 



momenta through the definitions 



ir (x,t) 

It- (x, t) 

n+(x,t) 



1 d 3 k 

~T d 3 k 

2u„ (2vr) 3 

~T~ d 3 k 

2u n (2vr) 3 
"T~ d 3 k 
(2vr) 3 



\a^e -- + a k e --} 



■{dfce*s + 4e~* 2 } 
(8) 



where 



a 6 a 6 a 6 a 2 

The quantum Hamiltonian was derived using 2 t± 

a Fourier mode decomposition of the fields and ^6 ~~ i m T. + 



4 



and we also define 



a 



6 



CO. 



7T± 



and 



= A[(< $ > 2 -v 2 ) + 2v 2 ] 
= A[(< $ > 2 -w 2 ) + 2v + t;_] 
= A[(< $ > 2 -V) + 2a 2 ] (10) 



For the rest of the section, for brevity, we drop 
the k,t dependence of the fi's and cu's. The 
Hamiltonian that describes the quantum evolu- 
tion of a DCC where the corresponding classi- 
cal condensate can be in any direction in isospin 
space, is: 



H H neu t ra i -\- Hcharged ~t~ Umixed 



(12) 



2v + v^+v 2 3 +a 2 = fl 



(11) where, 
l 



neutral 



d 3 k 1 u v , t t 
-{-j{a ] k a k + a k al) 



+ 



(2tt) 3 2 a 3 

2^ (^1 ~ 1 )(4 a fc + a fc°I + a -fe°fc + a -k a l) + -^1(44 + 44) 



+ 2^ ( !f " 1)(d * d * + 4 ^ fc + d " fcdfc + d ^ k)} 



(13) 



charged 



{-r(6fc6fc + Cfc4) + — s ( 



Co' 

(2tt) 3 L a 3 



2a 3 v 
I 



l)(6l6 fe + c4 + 6_ feCjfc + C t _ jfc 4)} (14) 



if, 



mixed 



I 



d 3 k \a 3 f 2 cos 2 {p)sin 2 {6) 



+4&* + c fc c_ fc + 4cL fc + c fc 4 + 4^ + 4 &t -fc) 

Aa 3 / 2 cos(p)sm(p)sm 2 (#) 

(b k a_ k + 6 fc a[ + 4afc + c k a_ k + c\.<J_ k + c k a[ + b[a k + 44 fc ) 

Aa 3 f 2 sin(p)sin(6)cos(6) / t t , t t \ 

(d fc a_ fc + 4< + 4 a * + d W-k) 

Aa 3 f 2 cos(p)sin(6)cos(8) 



M-fc + 6*4 + 44 + cfcGL fc + 44 fe + °kd k + 44 + 44 J} 



(15) 



5 



In order to show the dynamical origin of 
squeezed isospin states we consider the simplest 
case 9 = 0. In subsequent communications the 
other cases will be considered. 9 = implies the 
symmetry breaking takes place in the £ direction 
and the mixed term Hamiltonian , H mixe( i van- 



ishes. The general problem of evolution of the 
quantum state of the DCC with arbitrary orien- 
tation in isospin space will be considered later 
24l | . The total Hamiltonian, H for this special 



case reduces to 



f ^ — f k^TT / t t \ / \ 

+ 0(44 + 44) (17) 

+ ^3 ~ 1)(44 + 44 + d- k d k + dLA)} 

+ {-f (b{b k + c k c{) + - l){b{b k + c fe cl + b. k c k + cL fc 6[)} 



We notice that H has the form of a decoupled 
Hamiltonian . This is easy to understand from 
the ST) (4) parent. The SO (4) vector has been 
decomposed into four fields :7r±,7ro and £ being 
respectively the charged pions, the neutral pi- 
ons and the sigma fields. H has the character- 



istic quadratic su(l,l) structure that is associ- 
ated with Hamiltonians that are diagonalised by 
a squeezing transformation , which was explicitly 
done in 1. 

The parameters of the squeezing transforma- 
tion in the pion sector are 



I* = Cosh(r K ) = Jl[(9* +£.)+!] 
V 2 uj n il n 

and in the a sector are: 

I 



p = Cosh(r-z) = 




r and are the squeezing parameters, and p? — 
v 2 = 1 and p 2 — a 2 = 1. 



, = an*.) = ) /i[(£ + £)-!]. (18) 
I 



a = Sinh(r s ) = Jhp* + _ 1], (19) 

I 

Since the sigma field decouples in this par- 
ticular Hamiltonian, it can be analyzed inde- 



pendently of the pion fields. The diagonalized Hamiltonian can be written as: 
l 
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H 



d 3 k 1 rn ,. , t . 1, 



(ClC, + B\B t + !)} + a^D\D t + -)}. (20) 



where 



A k (t, r) = /i(r, t)a k + u(r, t)al k (21) 
D k (t, r) = p(r, t)d k + a(r, t)d[_ k 
C k (t, r) = /i(r, t)c k + u(r, t)b ] _ k 
B k (t, r) = fi(r, t)c_ fc + v(r, t)b\ 

Considering only the pion sector, the diago- 
nalized Hamiltonian H can be converted into a 
Hamiltonian in terms of quantum fields corre- 
sponding to the operators A,B,C and their ad- 
joints to obtain a purely quadratic Hamiltonian. 



H(t) 



(2vr) 3 ^ 2 vv a 

v ' i=A,B,C 



(22) 

The Schroedinger equation for each momen- 
tum mode is simply: 

H (k,t)^(k,t)=tj t ij(k,t). (23) 

If we use the n-representation (co-ordinate space 
representation) for tp(k,t), then, the su (1,1) 
symmetry of the Hamiltonian tells us that the 
solution for ip(k,t) is just a Gaussian. The 
equation satisfied by the wave functions for each 
mode are then given by: 



The expectation values of the number operator 
for the neutral pions for each momentum k is 
given by: 

< il) k (t)\a\a k \^ k (t) > = Sinh 2 (r) 

= < V> fc |4(*)4fc(*)W26) 
An identical expression holds for the charged 
scalar fields. 

The non-equilibrium transition is carried out 
as described earlier by making < $ > a function 
of time, where < $ > changes from to f n either 
instantaneously (quench) or adiabatically. In ref 
[i^, we have shown how the squeezing parameter 
r k is related to the competing effects of the ex- 
pansion rate of the plasma and the rolling down 
time of the system from a state of restored sym- 
metry to that of broken symmetry. We found 
that for the quenched limit (fast expansion) the 
low momentum modes are enhanced due to the 
squeezing parameter being very large, whereas 
for the adiabatic limit (slow expansion) no such 
enhancement occurs and the squeezing parame- 
ter is small. This enhancement corresponds to 
DCC formation. 



3a 



i/) A (k, t) + —i) A + (^Y(k, t)i> A (k, t) = 0. (24) 
a a 6 

(similar equations hold for fields B and C) where 



' a 



,, 2 (M) 



+ A(< $ > 2 (t) -v 2 ). (25) 



III. THE SU(1,1) DYNAMICAL SYMMETRY 
AND PION MULTIPLICITY DISTRIBUTIONS 



The su(l,l) symmetry of this Hamiltonian was 
studied in ref [23] . In terms of su(l, 1) generators 
the Hamiltonian can be written as 
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H 



d 3 k 1 
(2^)3^ 



2Q n ((fi 2 + v 2 )N + fw(V + X? f )) + 



d 3 k 1 
(2^)3^ 



(27) 

2fi s ((p 2 + a 2 )iV E + (7p(P E + £>*)). 



Where, satisfy a su(l, 1) algebra 

£> = a k a_ k + 6 fc c_ fc + c fc 6_ fc = K{ + + 
pt = a ^ a t + c t_ fc6 t + 6 t_ feC t = ^+ + x + + K + 

N = \i a \ a k + alk a -k + b[b k + bl k b- k + c\c k + c!_ fc c„ fc + 3} = K[ + K% + iT 3 ° 

(28) 



We also define: 



[N, V] = -V; [N, V ] ] = £>t; [V\ V] 



-2N 



(29) 



E 

Ny, 



dkd-k 
dl k 4 



\{d\d k + d ] _ k d_ k + 1} 



(30) 



which also satisfy an su(l,l) algebra. 

Thus, the Hamiltonian is linear in the genera- Shrodinger equation 
tors of the Lie Algebra of the SU(1, 1) group and 



following [27] the solution of the time dependent 



.#(£) 
dt 



d 3 k 1 
(2vr) 3 ^ 



2^((/i 2 + u 2 )N + ftu{V + 2^)^(0) 



(31) 



for each mode k , we take tanh(r k ) for both 
charged and neutral modes. Then we have 
\ip(t) >= Hk\if)(k,t) > and with discrete k we 
have for each mode k: 



i , = H k ip k [t) 



(32) 



, where 

H k = \2Q n (k,t)((fi 2 + u 2 )N + fiu(V + V^)). 
ar 

(33) 

The solution ip k (t) is given by the coherent state 
relevant to the discrete series representation of 
the SU(1,1) group given by: 



(34) 
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where, 

\a k (t) >= e <xk^ e nN kea >(k)v {l |^, (0) > _ (35) 

Here,afc = tanh(r(k)) and rjk = 
2ln(Cosh{r(k))) = -ln(l - \a k \ 2 ),a'(k) = -a%. 

r(k) is related to the frequencies fl n (k,t) and 
u v (k) by 

Tanh(2r k ) = f ■ - (36) 



We denote by no, n + and n_ the number of 
produced 7r s,7r^_s and n'_s. At non-zero mo- 
mentum, we have two such sets of number states, 
one for forward momentum (k) and one for back- 
ward momentum (-k). 

The states contributing to EHl are only the 
states which have the same number of modes i.e 



Where D, n (k,t — ► oo) = u n (k). Thus in the 
evolution of the condensate , it is the frequency 
changes which bring about squeezing [29t | . 



From the definitions of T>,T>< and N, For non- 
zero momenta, each of these exponentials has the 
form of a two-mode squeezed state which, for two 
generic modes a and b, is defined through [25]]: 



(37) 



(39) 



n ak = n a _ k , m bk = m c _ k and p Ck = p b _ k . For 
a complete statistical distribution of the pions 
we find the p-ordered characteristic function of 
the pion number distributions, this will allow us 
to construct all the correlation functions. The 
characteristic function of the neutral pion fields 
is 



o f„ \ „— at&t e (*<W tanh(r) „—(a^a-\-bb^)ln(cosh tanh(r) 

I 

with the two-mode squeezed state being given Using the above we get 
by: 

Kb >= S ab (a)\O a ; b > . (38) 
I 



|^(M) >= (sech(r k )) 3 a n k +m+p \n ak , n a _ k > \m bk ,m c _ k > \p Ck ,Pb_ k > 

n,m,p=Q 



Ci{£,V,p) =< afc ;0 a _JSt fcia _ fc (r)e« a l^^)e^-^ a -)^, a _ fc (r)|0 afc ;0 a _ fe > ePM'+M'AO ( 40 ) 

I 



field through the formula 

From this function it is relatively easy to con- 
struct all the "p-th ordered moments" of the pion 
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<* «fc a -k a -k a k > P~ 



(41) 



r 



. Similarly for the charged pions the Character- 
istic function is 



P(l«| 2 +l'?| 2 /2) 



and the correlation function is 



(42) 



< &V c m b n > - (—)i(—) l '~ S ^'~ S 



—r(—rc 2 (^, P )\^ 



r 



(43) 



Thus for example the forward backward correla- 
tion functions for neutral pions at non zero mo- 
menta are given by 

< a k a^ k >= sinh(r)cosh(r) =< a\,a}_ k > (44) 

Thus , there is an entanglement of the forward 
and backward neutral pions . 

< bkC-k >= sinh(r)cosh(r) =< b k c_ k > (45) 



< c k b_ k >= sinh(r)cosh(r) =< tf_ k c k > (46) 

This entanglement provides us with correlations 
between the pions with forward and backward 



momenta for both the charged and neutral sec- 
tor, the forward momentum positively charged 
pions with the backward momentum negatively 
charged pions and the backward momentum pos- 
itively charged pions and the forward momentum 
negatively charged pions. 

In the non-zero k limit however the distribu- 
tions of the neutral and charged pions are the 
same. This situation changes dramatically in the 
special case of soft pions k — ► 0. We will now 
show that in this limit the correlations of the 
neutral pions as well as the number distributions 
are different in the limit of large squeezing . To 
see this observe that in this limit the probability 
distribution splits up into the convenient form 



< n ,n + ,n.\^{t) k=0 > | 2 = | < n |e ro{a <> )2 - r °* a °|0 >< n+, n_ | e WcW* obC | > 



(47) 



where r = lim k >or(k) and 7q is the complex squeezing operator 

conjugate of r . Defining S(r ) as the one mode 
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which is just the product of squeezed distribu- 

, t , 2 t2 , tions for charged and neutral pions and only even 

S(m) = <nnle ro(a o ) r '* a o } l0> u t wl- 

u \' oj ^ oi l u ^ number of pions emerge. Writing n + = n_ = n c , 

= S n0i0 . (48) we get the distribution of charged particles to be 

S n+ , n_,o is then the two mode squeezing operator 



<n + ,n_\e^ btct - r *o bc ^\0 > 

^n-\-,n— ,0- 



(49) 



The neutral and charged pion distribution is: 

Pn ,n c =< S no fi > 2 < 5 ,t ^ + n _ )00 > 2 (50) 



P(n) 




r_0 =3 
r_0 =3 . 5 
r_0 =4 



50 100 150 200 



FIG. 1: Shows the variation of P„ c with n c for the ro = 3, 3.5 
and 4 



(tanh(r )) 2nc 
(cosh(r )) 2 



(51) 



While the distribution of neutral pions is 

n \(tanh(r )) m 



n ° ((^)!) 2 cosh(r o)2 n ° 



(52) 



P(n) 

.007 

0.006 

0.005 

0.004_ 

. 003^ 

.002 

0.001 



r_0 =3 
r_0=3.5 
r_0 =4 



*50 Too 150 200 n 



FIG. 2: Shows the variation of P no with no for the ro = 3, 3.5 
and 4 



Comparing the charged and neutral pion dis- responds to the total no. of pions < n >= 40 
tributions for a squeezing parameter which cor- 
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,(Fig. 3) and a squeezing parameter which cor- 
responds to < n >= 2400 (Fig.4)we see that 
there is a marked difference in the charged and 



neutral pion distributions for large values of the 
squeezing parameter for pions at zero relative 
momenta(soft pions). 



P(n) 




FIG. 3: Shows the comparison of P„ (n) and P„ c (n) for the 
< n >= 40 



In I we have found in the evolution of this 
Hamiltonian that in the quenched limit ( fast 
expansion) the low momentum modes are en- 
hanced significantly, whereas in the adiabatic 
(slow expansion) limit, no such enhancement oc- 
curs. This has been shown to be directly related 
to the value of the squeezing parameter, since 
for each mode < n >= Sinh 2 (r ). The differ- 
ence in the total charged and multiplicity distri- 
butions of pions produced at low momenta (soft 
pions) for large squeezing parameter is a direct 
mirror of this effect. Thus the measurable quan- 
tities like the charged and neutral total multi- 
plicity distributions at low momenta can be used 
to determine whether a quench or and adiabatic 
expansion occurred in a heavy ion collision and 
whether a disoriented chiral condensate formed 
or not. This serves as a very useful signal for 
non-equilibrium phase transitions in Heavy ion 
collisions. We now proceed to examine whether 
this dramatic effect survives theoretically the im- 
position of isospin conservation. 



IV. THE ISOSPIN SYMMETRY OF THE 
HAMILTONIAN 



We now examine the isospin structure of the 
pion Hamiltonian. Since isospin is conserved in 
strong interactions we need to extract states of 
fixed isospin from the squeezed number states 
which will be the eigenstates of the su(l,l) 
Hamiltonian constructed above. Thus far, most 
studies have constructed such isospin states at 
zero momentum. However, for an event by event 
analysis we need to construct these states at 
non-zero momentum states which have not been 
considered before. We will do so by defining 
the isospin operators constructed from the Fock 
state basis of pions at non zero momentum by 
constructing the vectors 



a(k) 



bk+Ck 
i(bk-Ck) 



\ 




(53) 



a(k) 



( fc-fc+c-fc \ 
i(b_fc-c_fc) 

a-k 



CLi( — k) 

a 2 (-k) 
a 3 (-k) 



(54) 



From the above ,two sets of isospin operators can 
be constructed 
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P( 
. 035 

. 03 
. 025 

. 02 
.015 

.01 
.005 



n) 



neutral 
charged 



FIG. 4: Shows the comparison of P no (n) and P na {n) for the < n >= 2400 



Ii(k) = ieijiaj(k)ai(k)^ 
Ii(-k) = ieijicijt-fyaii-ky . 



(55) 



su(2)+su(2)=so(4)defined by algebra the gener- 
ators 



Mi = Ii(k)+Ii(-k) 
N t = Ii(k) - Ii(—k), 



Generalizing the work of pj| , we find that the 
two operators 7j(/c)and h(—k) taken together de- 
fine a direct sum algebra of two su(2) algebras which satisfy the following relations 



(56) 



[M h Nj] = ieijM; [N h N,] = zey,M,; [M h M 3 ] = ie^Mi (57) 

I 

characteristic of an so (4) algebra. The Casimir which satisfy the nidified su(l,l) algebra: 
operators of this algebra are 



C\ = M 2 + N 2 = I(k) 2 + I(-k) 2 (58) 
C 2 = Mm N = I{k) 2 -I{-k) 2 

From the vectors a(k) and a(— k) we can form 
bilinear operators 



A(k) = a(k) ■ a{k) 
A ] {k) = otffc) • otffc) 
N(k) = a^{k) ■ a(k) 



[N(k),A(k)} = -A(k); 

[N(k),A(ky] = A(ky ; 

[A\k),A(k)} = -4N(k)-6 



(60) 



It is pertinent to note here that the relation- 
ship between the su(l, 1) algebra defined earlier 
(59) and the so (4) algebra is 



D=a(k)a(-k) = a 1 (—k)a\k) N = -(cft(k)a(k) + a\—k)a(—k) + 3) (61) 



Since the two algebras are related , we can also Since [H, M 3 ] = [H, N 3 ] = and 
relate the eigenstates of H constructed through [H, Ci]i=i,2 = 0, the eigenstates of M 3 , A^ 3 and 
these two algebras. Cj are defined by a generalization of the spheri- 
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cal harmonics Y^ l (a(ky)Y^ n (a^ — k^^O > which 
are also eigenstates of H. 



In terms of the Isospin operators the pion 
Hamiltonian can be written as: 



Hi, 



^-\2^ k {^ + v 2 )\{a\k)a{k) + a\-k)a{-k) + 3) + 2f2 7rk piv(a(k)a(-k) + a){-k)a\k)) 
[Ity)^ gt 2 



Since the pion Hamiltonian can be written both 
in terms of su(l,l) and isospin generators the 
eigenstates of the Hamiltonian have a Janus 
faced nature and can be both su(l,l) eigenstates 
and su(2)xsu(2) eigenstates. This allows the cal- 
culation of multiplicity distributions of pions in 
the number state basis as well as in the isospin 
basis as we shall see in the next section. 



ISOSPIN SQUEEZED STATES FOR k 
PIONS 



First, we construct the isospin squeezed states 
for the k — > case. For this we need to relate 
the states 



Fo >- 



; ro(X>t-X>) |,0 Q (O) >= e 7 -o(aJ 2 +4bJ+&o c o- a o+ 6 oco+co6o) ^(Q) 



> 



(63) 



r 



to the states of definite isospin. Note that we 
have added a subscript to the operators to in- 
dicate the zero momentum limit. 

For this purpose we note that for A; > 

there exists only one isospin algebra constructed 
from the operator 



% = 



/ &0+CQ \ 

i(b -c () ) 

\ a J 



This algebra is given in terms of 

I'iO — ^ijl a jO a lo 

The squared Isospin vector is 

I 2 = N 2 + N + AAo 



(64) 



(65) 



(66) 



For further discussion, now, in the k > 

case we drop the subscript . The definitions 
are 



a 1 • a 1 



A* 

A = a • a 

We also have the relations 



[N,A] 
[N,A*] 





-2A 
2A ] 



(67) 



(68) 



Yl(tf)\0> 



E(l-m)/2 
n=0 



;-2)- m v/(2/ + l)(/-m)!(/ + m)! 

om/2-ti 

{a ) l - m - 2n {a + ) n+m {a-) n \0 > 



(I — m — 2n)\n\(n + m)\ 



(69) 
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are the eigenstates of J 2 , I3 and N whose action 
on them is given by: 

\l,m,l + 2n> = N lin (A*) n Y^(rt)\0>(70) 

I 2 \l,m,l + 2n> = l(l + l)\l,m,l + 2n> 

I 3 \l,m,l + 2n> = m\l,m,l + 2n> 

N\l,m,l + 2n> = (I + 2n)\l, m, I + 2n > 



and the relationship between the number states 
\n , n + , n_ > and |Z, m, I + 2n > is given by 



n+(Z-H)/2 J, _ 1 1 

\l,m,l + 2n>= Yl cV' m ' n \p)\l + 2n-\m\-2p m+m + p, m m + p > (71) 



l—m 



i,™,n _ ( _ 1)m 2—+p (l + n) (1 + 2/)! (l + m)\n\ (I - m + 2n - 2p)\ (m + p)\ 
~ { } \l (/-m)!m! 2 (1 + 21 + 2n)\ (n - p)\ 2 p\ 

p (-iy 

j^(p-j)\ (n + j-p)\ (l-2j)\j\j\ 



r 



(72) 



The inverse relation is 

l 



n + +n_+n 

\n , n+, n_ >= c ( '' n+ ~ n ~ | ^++"j | + ™ ()) |Z, n + - n_, n + + n_ + n > (73) 

* — ' n + +n- =t-5 



For strong interactions involving charged pions only values of c^ m ' n contributing to physically 
n + — n_ = by charge conservation , hence the measurable quantities will be 

l 

i,o,n o(n ++ n- + fo > + OKgj + gKj + 2n - 2(n + + n_))!(n + + n_)!n!(n + + n_)! 1 

(»++»-) 1 (2n + 2/ + l)! J 



y («++«-) ( -*-) 

^ j " (2 2 ^')j! (-j + (n+ + n_))! ( 2./ • /)! (j - (n + + n_) + n)! 

r 



In considering the probability distribution of state of definite isospin has to be projected from 
pions with isospin conservation the squeezed 
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the number state distribution . Where we have assume that m > 0. 

n+(l-\m\)/2 

< n ,n + ,n^\S\l,m,2n >= c^' m ' n ^<n Q ,n + ,n^\S\l+2n—m—2p,m+p,+p> 

p= o The general expression for the one mode and 

(75) two mode matrix elements are 



q onemode(r) ,_ , s ^-m-2 P +n I + 2n - m - 2p\n \ tanhjr) ^ l+2u-m-2 P +n n 



mi n\ 



riQ l + 2n — m — 2p i 
2 ' 2 ' 



ur 9 ys T vJ (76) 

z_ (2A)!( + ~ m ~ p - A)!(^ - A)! 

I 



n,m, even. 



cTwo-mode _ / i N2n++m+p / j~ |7~ j ^l/^l (tanh(r)) n++m+P 
d n+ ,n-,m+p,p ~ { L ) ^ + ^ + V '^ V ' \cOsh(r)) n -- n ++ l 

( I U 

V A!(n+-A)!(m+p-A)!(n_ -n + + A)! 1 J 

A 



Using these identities we have 



n+(J-M)/2 

< n ,n_,n+|S(r)|Z,m,2n >= £ ^5^^^(r)52X^( r ) ( 75 

p=0 



For the special case of pion emission in strong To show clearly the effect of isospin conser- 
interactions , the charge is conserved, hence m = vation , we focus on total isospin / = as an 
n + — n_ = . example. In this case, 



< 0\S\n n n >- t^MT^l j ^° + ^ ! )h^( {2no)l ) (79) 
1 1 °' C (cosh(r))y* ( (2(n + n c ) + l)!) ) 2 ( K)! } ' (?J) 



r 



with n + + n_ = n c and n = tiq + n c . 



Thus, the probability of finding no neutral pi- 
ons and n c charged pions in the state of 1= 
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is This gives the neutral pion distribution P no 

(tanh(r)Y n ^ ((n Q + n c ))\ (2n )\ 

(cosh(r)Y l (2(n + n c ) + l)!) (n !) 2 

(80) as 

l 



2 -l-2n e tanh(r) 2 FC , sc ^ r) 2 + n ) _ F (l + n> tanh ( r )2)) ( sinh (r) 2 ) " (tanh(r) 2 ) ~ 2 " 

Pn ° ~ (cosh(r) 2 ) n n! 2 r(| + n) 

, (Ml 

In terms of the average number < n > of neutral 
pions 

l 



<ng> / \ i— no / \ 

2 -l-2n 0eTO ^ < no > -l + n ^ (2 no) , + n()) _ r( 1 + _^_ ) j 



P ™°" (l+<no>rno! 2 r(| + n ) 



r 



1M 



The coresponding distribution of charged pi- 
ons is given by 



2\ n 



\/7r li 7 ! ( 2 , 1 + n, tanh(r) ) Sech(r) (sinh(r) ) 



2(cosh(r) 2 ) n r(| + n) 1 ' 



r 



In terms of < n c > we have 



P - (<n>) " 2r(| + n c ) } (M) 

I 
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P( 

.002 




FIG. 5: Shows the variation P(n) with n with isospin I — 
for NEUTRAL PIONS for values of < n >= 500, 1500, 2500. 




FIG. 6: Shows the variation P(n) with n for CHARGED PI- 
ONS for values of < n >= 500, 1500,2500. 



To compare neutral and charged pion distri- 
butions for large and small squeezing parame- 
ters the following figures are presented. Fig. [6] 
shows the variation of P nc with n c for various 
values of < n >. Comparing the plots of the 
difference between charged and neutral pions dis- 
tributions without isospin (fig. 4) and those with 
isospin conservation (fig. 7), we find that the dra- 
matic difference in the two distributions is re- 
duced by the imposition of isospin conservation 



P(n) 




1 2 3 4 5 




12 3 4 



FIG. 8: Shows the variation P(n) with n for CHARGED PI- 
ONS (dashed line) and neutral pions (solid line) for values of 
ro = 5 

for 1 = 0. However, still the charged pion dis- 
tributions get broader with increasing squeezing 
and hence the distributions can still be used to 
determine whether a phase transition has oc- 
curred in Heavy ion collisions even when isospin 
conservation for soft pions is taken into account. 



FIG. 7: Shows the variation P(n) with n for CHARGED PI- 
ONS (dashed line) and neutral pions (solid line) for values of 
r = 3.5 
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VI. ISOSPIN SQUEEZED STATES FOR PIONS 
WITH ARBITRARY MOMENTUM 

In this section we calculate the multiplicity dis- 
tributions of squeezed pion states at non-zero 
momenta with the imposition of Isospin invari- 
ance. For arbitrary k we have to proceed slightly 
differently from above. We have seen that the 



states 

| a >= jVe oDt |0 > (85) 

represent the dynamical wavefunctions of the 
pions at non zero momenta where, D\k) = 
acft(—k)cft(k). Using the identitity: 



oo +1 

e k - = J2Mk 2 x 2 )J2 Y i m (kTYL(x), (86) 
1=0 -i 

I 



where 0; is the spherical Bessel Function <pi{x) = We write 

n + l)\x n 



and 



^H2'E (2 :;;% !nr w 

n=0 v ' 
I 



oo +1 

\a >= Ne °*(-k)*W\o >k | >_ fe =JV^^(a 2 #)U(-^) ]T ^(vM^) V^vM*^)!" >k \0 >- k 

1=0 m=-l 

m 



V 



To get the isospin decomposition we use the 
fact that 



ll(a(fc)t)|0,0,0> = (89) 
(-2)-™ v/(2/ + l)(/-m)!(/ + m)! 

Eit m)/2 a 2 T~ n u z ^ Ko)'" m " 2n K + ) n+m K-) n |o,o,o> • 

(/ — m — 2ny.nl(n + my. 



I 

y / m (a(/c)'l')|0 > is an eigenstate of l\ and 
/3fc and correspondingly Y[ n (a(— fc)^)|0 > is an 

eigenstate of I{-k) 2 and J 3 _ fc . |/,m,Z + 2n(fc) >= AU^r^VW)^ > fc 

The number operator eigenstates and the 

(90) 

isospin eigenstates are related by the formula: where N^ n = ( ( - 2n+2;+1 - )!n! - ) ) 2 
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Thus, the squeezed state wave function for 
fixed isospin is given by 



-io > k io >_ fc = > k io > fc 

(911 



n!(l + 2n + 2/)! 

n=0 m=-Z V ; 



Because, Y l _ m = (-l) m Y l rn (e)* and Y l m {ce) 
c'yjj(e) we get 



oo oo J a 2n+« 

la ' 



>= N i-n J2 J2 Yl 2 l (n + l)\ ( -~ 1 ^ l, ~ m,l + 2n>{k) l / ' m ' / + 2n> M) • ( 92 ) 



n=0 1=0 m- 



We now have the necessary ammunition to cal- and using the results of section III, we write 

culate the multiplicity distribution: the total probability of finding a state of 2n = 

non n 2 n(k) + n(—k) pions of momentum k and -k of 

Pi, m ,n + =< n o,n+,n-\k < n , n_, n + \^ k \a > . total Isospin 1= 1 and h = m as 

(93) 

Defining 



< n ,n + ,n_\ k \l,-m,l+2n > (k) = c l p m ' n (k) (94) 
and 

< n ,n+,n_|_ fc |/, -m,l + 2n >(- k )= c l ^ n (-k), 

(95) 



^ [a) - [ Sech(r)^(n(k)+n(-k)+l)\ Cp W ° p 1 )} W 



For 1=0 and m=0 we have 



n (±k),n.(±k)n + (±k) ( v (ta^(r)) 2 H fc )+"(- fc )+0 ((np(fc) + np(-fc) + n e (fc) + n c (-fc)))!2"°W+"»(- fc )+^' 

°'°' n W l 5ec/i(r) 3 (n(A;)+n(-A;)+/)! 1 (n (fc) + n (-k) + n c (fc) + n c (-k)) + 1)! 

, (M 



In figure 9 we show the effect of isospin con- multiplicity distribution of pions and see that 
servation at non zero momentum on the total 
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it has a significant effect of narrowing the multi- there is no difference between the neutral and 
plicity distribution. However as in the case with- squeezed pion distributions, 
out isospin conservation at non zero momentum 




FIG. 9: Shows the variation P(n) with n for pions with isospin conservation (1=0) (dashed line) and without isospin conservation 
( solid line) for values of ro = 5 



VII. CONCLUSION 

To conclude we have given a dynamical model 
for the production of pions in heavy ion collisions 
going though a non-equilibrium phase transition, 
and have shown how the charged and neutral 
pion distributions both globally and on an event 
by event basis can be used as a signal to mea- 
sure the nature of the phase transitions. We have 
shown that the distributions of the charged and 
neutral zero momentum pions are conspicuously 
different when the system undergoes a quench as 
is thought to occur when a disoriented chiral con- 
densate forms in a heavy ion collision. We have 
shown that this dramatic difference disappears 
when non zero momentum number distributions 
of pions are observed. However forward backard 
correlations are enchanced both in the neutral 
and charged pion sectors when the system un- 
dergoes a quench , which is parametrised by a 
squeezing parameter. We have shown how to in- 



corporate isospin conservation on both the zero 
momentum and the non-zero momentum pion 
distributions and shown that the effect of the 
quench survives the imposition of this invariance 
on the pion multiplicity distributions. Thus we 
have presented a field theoretical dynamical de- 
scription of the signals that can differentiate be- 
tween the various ways in which a system can 
undergo a phase transition in heavy ion colli- 
sions. The distributions and correlations which 
we have derived in this paper should be useful 
in testing the data envisaged at the RHIC col- 
lider to discern not only whether a disoriented 
chiral condensate is formed by easily measurable 
pion distributions, but also the effect of conser- 
vation laws on the general distributions of pions 



emerging in the collision processs[3l|. It is our 



hope that this work will prove useful to experi- 
mentalists when examining the data and we look 
forward to its application. 
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